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1. Introduction 


Consider a Hamiltonian system of 2 n differential equations 


(Xh) 


Qi — + 


8H 

8pi 


Pi 


8H 

8qi 


A first integral is a function of the qi and pi which is constant along the solutions of (Xh)- 
The system is called (meromorphically) Liouville integrable (or completely integrable ) when 
it admits n (meromorphic) first integrals F\..... F n which are functionally independent 
(their differentials are linearly independent) and in involution (their Poisson brackets vanish 
or, equivalently, the associated Hamiltonian vector fields Xp i commute). We refer to the ref¬ 
erence books |AM781 ICB97[ IAud08| for more on this topic; see also Section 2 for definitions. 


The Ziglin-Morales-Ramis theory (see |MRR10l IAud08| for statements and appli¬ 
cations) provides mathematical tools to check when a system is non-integrable. This 
is particularly useful as Hamiltonian systems generally come as parametrized families. 
The non-integrability criteria allow one to discard the vast majority of values of the 
parameters for which the system is not integrable. The principle is as follows. First, 
we find a particular solution T of the system ( Xh ) (generally from an invariant plane 
found from symmetries) and we compute variational equations (VE p ), i.e. systems of 
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linear differential equations governing a Taylor expansion of a solution of (Xh) along 
the particular solution T. The Liouville integrability of (Xu) induces integrability 
conditions on the variational equations (VE p ), which in turn imply properties of their 
monodromy or differential Galois groups. Technically, the Morales-Ramis-Simo theorem 
states that if ( Xh ) is integrable, then the Lie algebras of the differential Galois groups 
of all variational equations (VE p ) must be abelian (all these terms are defined in Section 2). 

The strength of this criterion is that it turns a geometric condition (integrability) into an 
algebraic one (abelianity of a Lie algebra), thus paving the way for possible computations. 
However, although there exist general algorithms to compute differential Galois groups 
of reducible systems such as the variational equations (VE p ) ( |Fenl5l lRet!4| or (vdH07| ). 
none of them are currently even close to being practical or implemented at this time. 
Furthermore, the size of the variational equations (VE p ) grows fast, so only a method 
which uses the structure of the system to make it simpler has a chance of being efficient. 
The main goal of the present paper is to explain how to use the structure of the system to 
make it simpler, which will allow us to check efficiently whether its Lie algebra is abelian 
or not. 

Over the past decade, several approaches have been devised to concretely apply this 
Morales-Ramis-Simo integrability criterion. 

For Hamiltonians of the form H = + F(g), where V is a potential in 

q, the first variational equation is often a direct sum of Lame equations of the form 
y"(x) = (n(n + l)p(x) + B ) y(x), where p denotes the Weierstrass function associated to 
an elliptic curve. In this case, Morales has elaborated a local criterion to find obstructions 
to integrability on higher variational equations via local computations (see Lemmas 11 and 
12 in |MRR01| Page 79, and Proposition 7, Page 81). Maciejewski, Przybylska and Du¬ 
val have elaborated techniques to handle variational equations for the case of Hamiltonians 
with potentials l |MP06l [DM091 [DM141 IDM15| L see also the works of Combot and coauthors 
|Coml31 IGKI21IBGSED14] , 

Another approach is to determine numerical trajectories and compute numerical mon- 
odromies around these. Although it is difficult to obtain rigorous proofs by these methods, 
they provide surprisingly precise information. They have been developed, for example, by 
Martinez and Simo [ MS09| . by Simon and Simo in the Atwood paper |PPR + 10] . by Simon 
in the more recent |Siml4al ISiml4bj and by Salnikov |Sall41 ISall3| . 

The general strategy for turning numerical evidence into rigorous proofs is to show 
that a certain commutator is non-zero. This in turn yields calculations of integrals and of 
residues, which can be achieved algorithmically due to their D-finiteness. This is used by 
Martinez and Simo in |MS09| and later systematized by Combot and coauthors, see e.g. 
[GKT 2 I ICombll IHGSFD1 4] , 

The approach that we develop in this paper follows previous work by two of the authors 
in [AMGW131 I AM Will, AMW12 ). We establish a reduction method. Consider the p-th 
variational equation (VE p ) : Y' = A(x)Y, where the coefficients of A(x) are in a differential 
field k. Given an invertible matrix P(x) (a gauge transformation matrix), performing the 
linear change of variable Z = P(x)Y produces an equivalent linear differential system for Z, 
denoted by Z' = P(x)[A(x)] Z. The principle of reduction methods is to look for a gauge 
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transformation P(x) such that the resulting system Z' = P(x)|44(x)] Z is “as simplified as 
possible”. 

Let G denote the differential Galois group of (VE p ) and g be the Lie algebra of G. 
Following traditional works of Kolchin and Kovacic, we will say that we have a reduced 
form when P(x)[A(x)] s g(k) (see Subsection I2.3.3[) . Despite the apparent technicality of 
this definition, the Kolchin-Kovacic theory shows why this is a desirable form. This is similar 
to the Lie-Vessiot-Guldberg theories of reduction of connections (see [BSMRlOl lBSMR.12] 
for the latter and their connections with the Kolchin-Kovacic theory of reduced forms). Our 
strategy in this paper is to compute such a reduction matrix P(x) efficiently. 

After this reduction process, the Lie algebra g is easily read and its abelianity (or not) 
is given in the process. Furthermore, if g is abelian, then this process will have prepared 
the system to allow an efficient reduction of the next variational equation. 


Our strategy can be summarized as follows. The p-th variational equation (VE p ) is a 
differential system of the form Y' = A{x)Y where A(x) has the form 


A(x ) 


A\{x) 

0 

S{x) 

A 2 (x) 


In the Morales-Ramis-Simo situation (see Subsection 1 2.5 [) , we may assume that the Ai{x) 
are in reduced form and that the Lie algebra of the differential Galois group of the block 
diagonal system 


y = ^diag^ 


with M diag = 


Ai(x) 

0 

0 

A 2 (x) 


has an abelian Lie algebra. We show (Theorem 13.31 in Subsection I3.2[) that the reduction 
matrix may be chosen of the form 


P[x) 


Id 

0 

Hih{x)Si 

Id 


where Id denotes the identity matrix, where the S) are easily found from S{x) and where 
the unknown functions fi(x) remain to be found. In Subsection 13.41 we show how standard 
linear algebra allows us to find these fi(x) as rational solutions of first order linear 
differential equations y' = A (x)y + C A( X ) where the c, are constant and where X(x) and 
the bi{x) are in a convenient field. 


Structure of the paper. In Section 2, we recall the necessary notions of Liouville 
integrability of Hamiltonian systems, differential Galois groups, reduced forms of linear 
differential systems and the Morales-Ramis-Simo integrability condition. This section 
contains only previously known material. In Section 3, we solve a problem that is 
interesting in its own right : given a block triangular differential system whose diagonal 
blocks are in reduced form and have an abelian Lie algebra, we give a practical procedure 
to put the system into reduced form (and hence compute its differential Galois group). In 
Section 4, we show how to reduce the Morales-Ramis-Simo condition to the latter problem 
and thereby provide an effective version of the Morales-Ramis-Simo integrability criterion. 
In Section 5, we demonstrate the efficiency of the method by computing the reduced 
form and the differential Galois groups of the first three variational equations on a four 
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dimensional example, originally considered in |CDMP10| . 
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2. The Morales-Ramis-Simo Integrability Condition 


2.1. Hamiltonian Systems and Liouville Integrability. Let (M, u>) be a complex 
analytic symplectic manifold of complex dimension 2n with n £ N*. Since M is locally 
isomorphic to an open connected domain U cz C 2ra , Darboux’s theorem allows us to choose 
a set of local coordinates (q , p) = (q\ ... q n , p\ .. . p n ) in which the symplectic form u 


is expressed as J : = 


o 

-Id„ 


Id„ 

0 


where Id T , 


denotes the identity matrix of size n. In 

these coordinates, given a function H £ C 2 {U ) : U —* C (the Hamiltonian), we define a 
Hamiltonian system over U <z C 2n as the differential equation given by the vector field 

Y , = TVH = v d JLA_ _ v d JLA_ 

H ’ “-j dpi dqi £-[ dqi dpi ’ 

corresponding to the Hamiltonian differential system 

(2.1) Qi = ^(q,p), Pi = -^(q,p), for * = 1... n. 

Consider a non-punctual integral curve T of (12.11) . A meromorphic function F : U —> C 
is called a meromorphic first integral of m along r if it is constant along integral curves 
in a neighborhood of T, or equivalently when Xh(F) = 0. Observe that the Hamiltonian 
is a first integral of ( 12 . 11 ) . as we clearly have Xh(H) = 0. 

The Poisson bracket {, } of two meromorphic functions f,g £ C 2 (U) is de¬ 
fined by {f, «} := (V f , JVo). In the Darboux coordinates, its expression is 

{f , g} = J] - pp. The Poisson bracket endows the set of first integrals with a 

i=l dpi UQi 

structure of Lie algebra. A function F is a first integral of (12.11) if and only if {F , H} = 0, 
i.e. H and F are in involution. Also, note that X^ F = [X F , Xh\, so the involution 
condition means that the associated Hamiltonian vector fields commute. 

A Hamiltonian system with n degrees of freedom is called Liouville integrahle by mero¬ 
morphic first integrals along the integral curve T if it possesses n first integrals (including 
the Hamiltonian) meromorphic over U which are functionally independent and in pairwise 
involution. 


2.2. Variational Equations. Among the various approaches to the study of meromorphic 
integrability of complex Hamiltonian systems, we choose a Ziglin-Morales-Ramis type of 
approach. Concretely, our starting points are the Morales-Ramis theorem |MRR01| and 
its generalization, the Morales-Ramis-Simo theorem |MRR,S071 IMRRIO] . These two results 
give necessary conditions for the meromorphic integrability of Hamiltonian systems. Here, 
we need to introduce the notion of variational equation of order p e N* along a non-punctual 
integral curve of m- 
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Let <h(.s, t) be the flow defined by the equation (12.11) . Given a non-punctual integral curve 
r of (12.11) and zo e T, we let cj){t) : = j t) denote a temporal parametrization of L. We 
define the p th variational equation (VE^) of (12.11) along T to be the differential equation 

satisfied by the := - ’ ^ ^ for j p. For instance, the first three variational equations 
are given by (see |MRRS07| . §3.4, Equation (14), Page 860): 


(VEj) : 


(VEj) : Cl = dpXni, 

i'2 = d^X H (£ 1 , £i) + d(jjX H t;2 
?3 = d^X H {i\ , 6 , 6 ) + 2d^X H (£i , £ 2 ) + d^XnC: 3. 


(VEj) : 

For p = 1, the first variational equation (VE^) is a linear differential equation 


£1 = A\^i where A\ := d^Xn = J ■ Hess^H) e sp(n, C(<fo(t))), 


where C(</>(t)) denotes the differential field generated by the coefficients of the parametriza¬ 
tion (f)(t). Higher order variational equations are not linear for p ^ 2. However, for every 
(VEJ), one can construct an equivalent linear differential system (LVE^) called the lin¬ 
earized p th variational equation (see [ MRRS07] . §3.4 and |Siml4bj ). Indeed, (VE^) is linear 
in £ p and polynomial in the for i < p\ however, the for i < p are solutions of the linear 
differential system (LVE^ *) so that polynomials in the also satisfy linear differential 
systems, obtained via symmetric powers and tensor constructions. See, for example, §3 of 
|AMCWl3] for practical details on these tensor constructions on differential systems. 

For example, (VE?) is linear in £2 and linear in the monomials of degree 2 in the £ 1 , i.e. in 
the solutions of the second symmetric power system Y' = sym 2 (Ai)y. Hence the system 
(LVE^) is lower block-triangular and its diagonal blocks are sym 2 (Ai) and A\. We obtain 
(see e.g. |MRR10l I AMW1 H ISiml4bl ICW 15] ) the following matrices A p for the first (LVEp: 


sym 2 (Atix)) 

0 

S 2 {x) 

Ai{x) 


sy^{A l {x)) 

0 

Sz{x) 

A 2 {x) 


A 3 (x) = 

V 

In general, the matrix of (LVE^) is of the form 

A p (x) = 


sym 3 (Ar(x)) 

0 

0 

S 3 , 2 (®) 

sym* (Ai(x)) 

0 

S 3 , i(*) 

S 2 (x) 

Ai(x) 


sym p (Ar(x)) 

0 

S p (x) 

A p - 1 (x) 


In |Siml4b] , §4.1, Simon provides explicit formulas for these linearized variational equations. 
In what follows, we will identify (VE(*) and (LVE^) and we will just speak of variational 
equations of order p. 

The matrix sym 4 (Ai(x)) has ( Tl +^ 1 ) rows and columns, so that (LVE^) is a first order linear 
differential system of d p : = ^f =1 = ( n+p ) — 1 equations. The size d p grows fairly 

fast (polynomially of degree n in p) and forbids the use of a generic algorithm to compute 
on (LVEp. For this reason, we elaborate a specific algorithm which takes advantage of 
the structure of (LVE^) so that the polynomial growth of the size will become a relatively 
minor concern. 
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2.3. Differential Galois Theory and Reduced Forms. We begin this subsection with 
elements of differential Galois theory. We refer to |PS03| or |CH1 11 ISin09| for details and 
proofs. 

2.3.1. The Base Field. Our base field will be k := C(</>(f)), the differential field generated 
by the coefficients of the parametrization (j>{t ) (and C is the field of constants, which is 
assumed to be algebraically closed). We need to make assumptions about k to elaborate 
our algorithms. First we assume that k is an effective field, i.e. that one can compute 
representatives of the four operations +,—, x,/ and one can effectively test whether two 
elements of k are equal, see e.g. [Sin nu. Secondly, we assume that, given any scalar linear 
differential equation L(y(x)) = 0 where 

L(y(x)) : = a n (x)y^ l ' > (x) + a n _i(x)y (n-1) (x) H- \-a 1 (x)y , (x) + a 0 (x)y(x), with ai(x) £ k, 

one can effectively compute a basis of its space of rational solutions, i.e. the solutions 
which are in the base field k. The standard example of such a held would be k = C(x ) 
with C = Q. Singer showed, in |Sin91| . Lemma 3.5, that if k is an elementary extension 
of C(x ) or if k is an algebraic extension of a purely transcendental Liouvillian extension 
of C(x), then k satisfies the above two conditions and hence suits our purposes. He also 
proved, see Theorem 4.1 in (Sin9lj . that an algebraic extension of k still satisfies our two 
assumptions, which will be useful, as reducing the first variational equation may induce 
algebraic extensions. 

2.3.2. Differential Galois Theory. Let us consider a linear differential system of the form 
Y'(x) = A(x)Y(x), with A(x) £ A4 n (k), that is a square matrix of size n e N* in 
coefficients in k. A Picard-Vessiot extension for Y'(x) = A(x)Y(x) is a differential held 
extension K |k, generated over k by the entries of a fundamental solution matrix and such 
that the held of constants of K is C. The Picard-Vessiot extension K exists and is unique 
up to differential held isomorphism. 

The differential Galois group G of Y'(x) = A[x)Y ( x ) is the group of held automorphisms 
of the Picard-Vessiot extension K that commute with the derivation and leave all elements 
of k invariant. Let U(x ) £ GL n (iL) be a fundamental solution matrix of Y'(x) = A(x)Y(x) 
with coefficients in K. For any <p £ G, (p(U(x)) is also a fundamental solution matrix, so 
there exists a constant matrix C^ £ GL n (C) such that ip(U(x)) = U(x).C lf . The map 
pu : ip i—> is an injective group morphism. An important fact is that G, identihed 

with Im pu, may be viewed as a linear algebraic subgroup of GL n (C). 

Two linear differential equations Y'(x) = A(x)Y(x) and Y'{x) = B(x)Y(x), with 
A(x),B(x) £ A4 n (k) are said to be equivalent over k (or gauge equivalent over k) when 
there exists P{x) £ GL n (k), called a gauge transformation matrix, such that 

B(x) = P(x) [4f(x)] := P(x)A(x)P~ 1 (x) + P , (x)P~ 1 (x). 

Note that in this case: 

Y\x) = A[x)Y[x) <=> [P(x)Y(x)] > = B(x)P(x)Y{x). 

Conversely, if there exist matrices A(x), B(x) £ A4 n (k) and P(x) £ GL n (k), such that we 
have Y'{x) = A{x)Y{x), Z'{x) = B{x)Z{x) and Z{x) = P{x)Y{x), then 

B(x) = P(x) [A(x)J . 
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The Lie algebra g of the linear algebraic group G cz GL n (C) is the tangent space to G 
at the identity. Equivalently, it is the set of matrices N such that Id n + eN satisfies the 
defining equations of the algebraic group G modulo e 2 . 

Part two of the following proposition is known as the Kolchin-Kovacic reduction theorem. 
A proof can be found in |PS03j . Proposition 1.31 and Corollary 1.32. See also |BSMR10| . 
Theorem 5.8. 

Proposition 2.1 (Kolchin-Kovacic reduction theorem). Let us consider the differential 
system Y'(x) = A{x)Y{x) with A{x) £ _A4 n (k). Let G be its differential Galois group and g 
be the Lie algebra of G. 

(1) Let H cz GL n (C) be a linear algebraic group and f) cz _A/f n (C) be its Lie algebra. If 
A{x) belongs to f](k) := 1) (xwk, then G is contained in a conjugate of H. 

(2) Assume that k is a C l -field\j\ and G is connected. Let H zd G be a connected linear 
algebraic group with Lie algebra f) such that A(x) 6 f](k). Then, there exists a gauge 
transformation P(x) £ H( k) such that P(x)[A(x)] £ g(k). 

2.3.3. Reduced Forms of Linear Differential Systems. Let A{x) £ Ad n (k), G be the differ¬ 
ential Galois group of Y'(x) = A(x)Y(x) and g its Lie algebra. 

We say that the system Y'(x) = A(x)Y(x) is in reduced form (or in Kolchin-Kovacic re¬ 
duced form) when A{x) £ g(k) = g®^-k. This section contains a quick survey on reduced 
forms and their practical use. 

Following |WN63| . a Wei-Norman decomposition of A{x) is a finite sum of the form 

A(x) = Yj a i(x)Mi , 

where Mi has coefficients in C and the afix) £ k form a basis of the C-vector space spanned 
by the entries of A{x). The Mi depend on the choice of afix) but the C-vector space 
generated by the M t is independent of the choice of the afix). 

Definition 2.2. Let Lie(M) c Ai n (C) denote the Lie algebra generated by the Mi. We 
define Lie a i g (A) cz A4 n (C), called the Lie algebra associated to A, as the algebraic envelope 
of the Lie algebra Lie(-A), i.e. as the smallest Lie algebra of a linear algebraic group which 
contains Lie(A). 

Let Lie(A;k) := Lie(A)(k) cz M n ( k) and Lie a i g (A;k) := Lie a i g (A)(k) cz M n ( k). We 
see that the system Y'(x) = A{x)Y{x) is in reduced form when Lie a i g (A;k) = g(k). 

These reduced forms have long been studied in the context of inverse problems in 
differential Galois theory (see |MS02| and references therein). Their use in direct problems 
is more recent. Blazquez and Morales use them in their studies of Lie-Vessiot systems in 
(BSMRIOi IBSMR12] . Their application to Morales-Ramis theory is initiated in |AMW12| 
where Aparicio-Monforte and Weil show how to put the first variational equation in 
reduced form. In |AMCW13j . the same authors with Compoint show that a system is 
in reduced form if and only if, for any tensor construction const(A(x)) on A{x), any 
rational or hyperexponential solution of Y' = const(M(x))K has constant coefficients. 


*A field k is a C 1 -field when every non-constant homogeneous polynomial P over k has a non-trivial 
zero provided that the number of its variables is more than its degree. For example, C(x) is a C 1 -field and 
any algebraic extension of a C 1 -field is a (T-field. 
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One can also find in |AMCW13| a complete procedure to put a linear differential sys¬ 
tem into reduced form when it is irreducible (or completely reducible). This does not 
really apply here as the variational equations are generally reducible (and not completely 
reducible) systems - and the reduction method of [AMCW13] is far from being efficient yet. 

The approach that we elaborate in this paper was initiated (incompletely) in |AMW11| . 
It is based on another criterion for reduced form, which is given in the following lemma. 

Lemma 2.3. Given A(x) £ _M n (k), let G be the differential Galois group ofY'(x) = 
A{x)Y{x) and g be its Lie algebra. Let H be a connected linear algebraic group whose Lie 
algebra f) satisfies ti = Lie a i g (A). Assume that G is connected. 

Then Y'(x) = A(x)Y(x) is in reduced form, i.e. G = H and Q = t), if and only if, for all 
gauge transformation matrices P(x) in H( k), we have f)(k) = Lie a i g (P[A]; k). 

Proof. Follows directly from the Kolchin-Kovacic reduction theorem, see Proposition 12.11 

□ 


2.4. The Morales-Ramis-Simo Integrability Criterion. We are now in position to 
state the Morales-Ramis-Simo integrability criterion. See |MRRS07] for a proof and Ij2] for 
the definitions. 

Theorem 2.4 (Morales-Ramis-Simo integrability criterion). Consider a Hamiltonian vector 
field Xu and a non-punctual integral curve T. For p £ N*, let G p be the differential Galois 
group of (VEp, the p th variational equation along T. Let g p be the Lie algebra of G p . 
Assume that the Hamiltonian vector field Xu is Liouville integrable by meromorphic first 
integrals along the integral curve T. Then, for all p £ N* ; g p is abelian. 

Of course, given p £ N*, computing the differential Galois group G p of such a big dif¬ 
ferential system would be an unrealistic task in practice unless we use the structure of the 
system to simplify the computations. We will establish a specific reduction method , i.e. 
compute a gauge transformation matrix P p (x ) such that P p (x)[A p (x)\ £ g p (k). After this 
reduction process, the Lie algebra g p is easily read and its abelianity (or not) is given in 
the process. Furthermore, if Q p is abelian, then this process will have prepared the system 
to allow an efficient reduction of the next variational equation. 


2.5. The Strategy for an effective Morales-Ramis-Simo Criterion. We refer to *12.21 
and *12.31 for the notations used in this subsection. Let us fix an integer p ^ 2. The matrix 
of the p th variational equation has the form 


A p (x) 


sym^M^x)) 

0 

S p (x) 

Ap-i(x) 


For each m £ {1,... ,p}, we let G m denote the differential Galois group of the m th variational 
equation Y'(x) = A m {x)Y{x) and g m its Lie algebra. For all m £ {1 ,... ,p— 1}, we assume 
that we know a gauge transformation matrix P m (x) such that P m (x)[A m (x)\ is in reduced 
form, i.e. Lie a i g (P m [A m ]) = g m , and we further assume that each g m is abelian. We let 
A m red{ x ) denote the obtained reduced form, that is A m}re d(x) := P m (x)[A m (x)\. 

Under these hypotheses, we will show in the next section how to put the p th variational 
equation A p [x) into reduced form in an efficient way. 
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Remark 2.5. Our assumptions imply that the first variational equation is in reduced form. 
This in turn implies that our base field k is no longer just C(fj) but may be an algebraic 
extension of the latter (see |AMCW13| ). In the sequel, our base field k is the algebraic 
extension of Cfifj which is needed to put the first variational equation into reduced form. 
Since an algebraic extension of a C 1 -field is a C 1 -field, we obtain that k is a C 1 -field provided 
that Cfcfi) is a C 1 -field. Consequently, we are allowed to use Proposition \2. II as soon as Cffij 
is a C 1 -field. From now on, we assume that k is a C 1 -field. 

Our assumptions also imply (see |AMCW13] . Lemma 32, Page 1513) that, for all 
me {1,... ,p — 1}, the differential Galois groups G m are connected. Moreover, both the 
groups G m and their Lie algebras Q m are abelian. 

Lemma 2.6. The group G p is connected. 

Proof. This is a direct application of |MRR10| . Lemma 10. □ 

As we can see in |AMCW13j . Lemma 14, Page 1508, 

Sym p (Pi(x))[sym p (Ai(x))] = sym p (Ai ire(i (a:)). 

Also, sym p (Ai )T . e d(x)) is a reduced form of sym p (Ai(x)). Indeed, this follows from 
|AMCW13| . Theorem 1, because any tensor construction on sym p (Ai(x)) is a construc¬ 
tion on A\{x). 

Consider the block-diagonal gauge transformation matrix 


Sym p (Pi(x)) 

0 

0 

P p -i(x ) 


Q{x) := 

Thanks to the above remarks (see also |AM10| . §4.5.2), we find that 
Q(x)[A p (x)] = 


sym p (Al,red(®)) 

0 

S(x) 

Ap—i ,red(pE) 


where S(x ) has entries in k, and the block-diagonal part of Q(x)[A p (x )] is in reduced form. 


Furthermore, Lie a i g 


sym p (Agred) 

0 

0 

l,red 


is abelian. 


3. Reduction of Linear Differential Systems with a Reduced Abelian 

Diagonal Part 

The previous subsection shows that finding a reduced form for the p th variational equation 
now amounts to finding a reduced form for 


A{x) := Q(x)[A p (x )] = 


sym p {A ltred (x)) 

0 

S(x) 

Ap— l J red(*^') 


e .Mn(k). 


The submatrices sym p {Ai red (x)) and A p _i tred (x) belong respectively to A4 ni (k) and 
A4 na (k), with 77,1 := and n 2 := ( n+ ^ 1 ) - L 


We have A(x) = A diag (x) + A sub (a;), where A diag (x) := ^ Sym (-di.rcvD)) 


sym p (Ai tred (x)) 

0 

0 

Ap— l 5 red(^') 


and A sub (x) := 


0 

0 

S{x) 

0 

is abelian. 


We have seen that Y\x ) = Adiag(®)^(®) is in reduced form 


to put the full system Y'(x) = A(x)Y(x) in reduced form. 
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3.1. The Diagonal and Off-Diagonal Subalgebras. We refer to £12.3.31 for the no¬ 
tations used in this subsection. Let £ M. n {C) be a basis of Lie a i g (Adiag) 

and let £ _A/l n (C) be a basis of Lie a i g (A su b). We define the vector space 

tl := Lie a i g (^diag) ® Lie a i g (Al su b). Note that Lie a i g (A) £ fp and 1) is the Lie algebra of 
a linear algebraic group. Let us sum up some elementary properties of f) in the two follow¬ 
ing lemmas: 


Lemma 3.1. Let us consider a matrix 


£ Lie a i g (j4 su b; k). 


iVi(x) 

0 

N 2j \(x) 

N 2 (x) 


e fj(k) 


0 


f ° 

0 

C\(x) 

o ) 

1 C 2 (x) 

0 


(1) For (i,j) £ {1; 2} 2 , 


(2) The matrix 


0 


( 0 

0 

Ci(x) 

0 ) 

V Cj(x) 

0 


= 0 . 


iVi(x) 


mix) 

° 

f ° 

0 

N 2 , i(x) 

N 2 (x) ) 


0 


and the 


N 2 ,i(x) 


N 2 {x) 


;ie a i g (ffsubi k) is an ideal in f)(k 


0 

0 

Ci(x) 

0 

l in f)(k 



belong to Lie alg (A sub ; k). 


and matrices 


Lie bracket 
Furthermore 


Proof. (1) A straightforward computation shows the first point of the lemma. 


(2) We have 


N^x) 

0 

N 2t i(x) 

N 2 (x) 


Ci(x) 


0 


N 2 {x)Ci(x) 


0 


£ f)(k) and 


n^x) 

° ^ 

f ° 

0 V 

( 0 

0 

N 2 ,i(x) 

N 2 (x) ) 

V mix) 

0 )\ 

V N 2 (x)Ci(x) - Ci(x)N 1 (x) 

0 


efl(k). 


We prove that they belong to Lie a i g (A su b; k) using that fact that the diagonal 
blocks of the two matrices are zero. The latter Lie bracket identity also shows that 
Lie a i g (A su b; k) is an ideal in f)(k). 

□ 


Lemma 3.2. For all B(x) £ Lie a i g (A su b; k), we have exp (B(x)) = Id n + B(x) and 
log(Id n + B(x)) = B(x). This induces two bijective maps which are inverses of each other 

exp . Lie a i g (A su bi k) > ^Id n + B(x'), B(x') £ Lie a j g (-A S ub> k) 1- 

B(x) >-> Id n + B(x) 

log ■ |ld n + B(x), B(x) £ Lie a i g (A su b, k) > Lie a i g (A su b, k) 

Id n + B(x) >-> B(x). 

Proof. Let B{x) £ Lie a i g (A su b; k). The equality exp(R(x)) = Id n + B(x) is a di¬ 
rect consequence of the first point of Lemma 13.11 The same argument shows that 
log(Id n + B{x)) = B(x). It follows directly that exp and log are bijective on the wished 
sets and inverses of each other. □ 


3.2. The Shape of the Reduction Matrix. We refer to £ 12.31 and £ 13. II for the notations 
and definitions used in this subsection. The aim of this subsection is to prove: 
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Theorem 3.3. There exists a gauge transformation 

P{x) 6 |ld n + B(x),B(x) £ Lie a i g (A sub ; k) j, 
such that Y'(x) = P(x)[A(x)]Y (x) is in reduced form. 

Let G be the differential Galois group of Y'{x) = A{x)Y{x). By construction, we have 
G = G p , where G p is the differential Galois group of the p-th variational equation Y'(x) = 
Ap(x)Y{x). Let H be the connected linear algebraic group with Lie algebra f). Before 
proving Theorem 13.31 we start with a key lemma. 


Lemma 3.4. There exists a unipotent gauge transformation P(x), of the form P(x) = 
£ H{ k), such that Y'[x) = P(x)\A(x)]Y(x) is in reduced form. 


Id ni 

0 

N(x) 

Id n2 


Proof. Let Ha be the connected linear algebraic group with Lie algebra Lie a i g (A; k). 
We have the inclusions G c Ha Q if. As G = G pi Lemma 12.61 shows that G is 
connected. So we may use the second point of Proposition 12.11 to obtain the exis¬ 


tence of Q(x) : = 


Di(x) 

0 

Sq(x) 

D 2 (x) 


£ Ha( k) such that the linear differential system 


Y'(x) = Q(x)[A(x)]Y(x) is in reduced form. Let R(x) := 


Df\x) 

0 

0 

Df\x) 


£iL(k) 


Id ni 

0 

d 2 \x)Sq(x) 

Id n2 


£ H{ k). Consequently, to prove the lemma, 
it is sufficient to prove that Y\x) = R(x)Q(x)[A(x)]Y(x) is in reduced form. We have to 


prove that Lie a i g ^Q[A];kJ = Lie a i g ^i?Q[A];kJ. Let H r q be the algebraic group whose 
Lie algebra is Lie a i g . Thanks to the first point of Proposition 12. 11 the group H R q 

contains G. Since Y'(x) = Q(x)[A(x)]Y(x) is in reduced form, G is an algebraic group 
whose Lie algebra is Lie a i g . This implies that Lie a i g kj c Lie a i g k 


Let Ji|k denote the Picard-Vessiot extension for the equation Y'(x) = A(x)Y(x) 

£ GL n (iL), with Ui(x) e GL 11; (K) be a fundamen- 


U^x) 

0 

U 2 ,l(x) 

U 2 (x) 


and let U(x) : = 
tal solution. The elements of G are of the form 
Gi £ 


G\ 

0 

G2, 1 

g 2 


GL n (C), with 


Id ni 

0 

G 2 ,i 

Idn 2 


GL n . (C). Let G sub be the subgroup of elements of G of the form ( ^ ni 

\ ^ 2,1 

A direct computation shows that G sub is a normal subgroup of G. Therefore, G ~ 
G sub x G/G su b- Due to [PS03| . Proposition 1.34, (2), Gviia. g := G/G sub is isomorphic 
to the differential Galois group of Y'(x) = Afc ag (x)Y (x). Let us write Q{x)[A(x)~\ =: 


Di(x)[sym p (A lired (x))] 

0 

—2,1 (*^) 

-^2 («e) l,red(*^)] 


in k. We use the relation G ~ G aub x Gdi ag and the fact that Y'(x) = Q(x)[A(x)]Y(x) is 
in reduced form to find that 


Lie. 


alg 


(Q[A]-, ] 


k) ~ Lie 


alg 


D\ [sym p {A]^ e< i)\ 

° 1 

(k) © Lie a i g 

( 0 

0 

0 

D 2 1 ,red\ j 

l —2,1 

0 
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A direct computation shows that 
(3.1) R(x)Q(x)[A(x)] = 

By construction, 


sym p (A^ red (x)) 

0 

Df\x)A 21 {x)D 1 (x) 

Ap— l 5 red(^) 


Lie a i g ( RQ[A]; k) c Lie, 


-alg 


sym p {Ai^ red ) 

0 1 

(k) ® Lie a i g 

( o 

0 

0 

Ap—i^red j 

V L> 2 ^2,1^1 

0 


Lie a i g 

Lie a i g 

(3.2) 


Since D\(x) and D 2 (x) are invertible matrices, Lie a i g ^ —-(k) and 

(k) have the same dimension. 


0 

0 

L*2 1 ^2,1-L ) 1 

0 


(Q[A];k) c Lie a i g (RQ[A]; k) we obtain that 


Lie a i g 


0 

0 

— 2,1 

0 


(k) = Lie ; 


alg 


0 

0 

v a 2)X d x 

0 


Due to the inclusion 


(lO- 


Using the facts that the systems Y'(x) = Adia.g(x)Y (x) and Y'(x) = Q(x)\_A(x)~\Y (x) are 
in reduced form and G ~ G su b x Gdiagj we find that 


Lie. 


alg 


sym p (Ai red ) 

0 1 

(k) = Lie alg 

( T>i[sym p (Ai, red )\ 

0 

0 

l,red J 

K o 

D 2 [-*4p— 1 ,red\ 


(k). 


Combined with (13. 2p . this proves that Lie a i g \RQ(A\] kj cz Lie a i g ^Q[A]; k^. Since 


we 


have an inclusion Lie a i g kj c Lie a i g kj, we obtain the equality 

^Q[A];kV In other words, Y'(x ) = R(x)Q(x)[A(x)]Y(x) is 

□ 


Lieaig (^Q[A];k) = Lie, 
in reduced form 


Proof of Theorem \3.A It follows from Lemma 13.41 that a reduction matrix can always be 


chosen of the form P(x) = 


Id ni 

0 

N(x) 

Id n2 


£ H( k), where N(x) £ ,M n2 . ni (k). By a 


straightforward computation, we find log(P(x)) = 


0 


N(x) 


0 


0 


£ fl(k). But with the 


same reasoning as in the proof of Lemma 13.11 we obtain that log (P(x)) £ Lie a i g (A su b;k). 
This concludes the proof of Theorem 13.31 □ 


The following corollary will be crucial for the reduction procedure of 

Corollary 3.5. Assume that, for all gauge transformations of the form 
P(x) £ jld n + B(x), B(x) £ Lie a i g (A su b;k) j, we have Lie(A;k) = Lie(P[A];k). Then, 
Y'(x ) = A(x)Y(x) is in reduced form.. 

Proof. Theorem 13.31 provides a B{x) £ Lie a i g (A su b; k) and P(x) = Id n + B(x) such that 
the system Y'(x) = P(x)[A(x)]Y(x) is in reduced form. In virtue of the hypothesis, 
Lie(A;k) = Lie(P[A];k). This implies that Lie a i g (A;k) = g(k), where g is the Lie algebra 
of the differential Galois group G of Y'(x) = A(x)Y(x). This proves that Y'(x) = A(x)Y(x) 
is in reduced form. □ 
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3.3. The Adjoint Action. We refer to H2.3I and iJ3.ll for the notations and definitions 
used in this subsection. In 1)3.21 we have proved the existence of a gauge transformation 
matrix P(x) e jld n + B(x), B{x) e Lie a i g (A sub ; k) such that Y'(x) = P{x)\_A{x)~\Y(x) is 
in reduced form. Let B\,..., B a 6 A4 n (C) denote a basis of Lie a i g (A sub ). 

a 

Proposition 3.6. If P(x) := Id n + ^ fi(x)Bi, with fi(x) e k and Bi e Lie a i g (A sub ), then 

1 = 1 

<j a 

P{x)[A{x)] = A(x ) + Y fi(x)[Bi, A diag (z)] - Yj 

i=1 i=L 

Proof. Due to the first point of Lemma 13.11 we have the equalities 

a a 

P~Hx) = Id n - Y fi( x ) B i and P(x)A(x) = A(x) + Y fi{x) B iA d iag(x). As 

i= 1 i=1 

A(x) = Adiag(ic) + ArabC®)) we use Lemma 13.11 and find that 


P{x)A(x)P 1 (x) = ( A diag 0) + A sub (x) + Y fj( X ) B j A diagO) I ( Id n “ Y fk( X ) B k 


3 = 1 


k= 1 


At^x) + fj{^x)Bj A d i ag (x) "y ' A d j ag (x)5fc 


3 = 1 


k= 1 


A(x) + Y fi{ x )[ B i, A diag(z)]- 


i =1 


Similarly, we have 


P\x)P-\x) = ( 2 fl(x)Bij |M n - YfjWBj ) = S 


\i=l 


i= 1 


This yields the desired result. 


□ 


We have seen in Lemma EO that Lie a i g (A sub ; k) is an ideal in fj (k). In particular, for 
all B(x) e Lie a i g (A sub ; k), the bracket [B(x), A diag (x)] is in Lie a i g (A sub ; k). This implies 
that the k-linear map T := [•, A d i ag (a:)], which is the adjoint action of Lie a i g (A d ; ag ;k) on 
Lie a i g (A sub ; k), is well defined: 


'L . Lie a i g (A sub , k) > Lie a j g (A sub , k) 

B(x) i —> [B(x),A diag (x)]. 

The following lemma will be necessary in 1)3.41 Note that the proof of the lemma gives 
a complete description of a finite set containing the eigenvalues of T. 

Lemma 3.7. The eigenvalues of the linear map 'll belong to k. 

Furthermore, there exists a basis of constant matrices, such that the matrix of the linear 
map T in this basis is block-diagonal, with blocks that are upper-triangular matrices with 
only one eigenvalue. 

Proof. Let M±, ..., Ms £ A4 n (C) be a basis of Lie a i g (A d j ag ), which is abelian. We may write 
4l<liag(.r) S gi(x)Mi with gi(x) e k. Let VH, := [•, M,\ denote the adjoint action of Mi on 

i= 1 
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Lie a i g (A su b). As the matrices Mj commute pairwise, the Jacobi identity on Lie brackets 
implies that the 'Ll also commute pairwise. The T,; have coefficients in the algebraically 
closed field C and commute pairwise, and therefore they are simultaneously triangularizable 
in a basis (Cj) of Lie a i g (A su b ). By construction, the Cj are constant matrices. Each Cj 
lies in a characteristic space of 47 associated with an eigenvalue \ t . 3 . We define Xj(x) : = 

2i=i 9i( x )Xij- As T = Xii=i 9i( x )^ii we see that the A j(x) £ k are the eigenvalues of *k 
and that the matrix of T is triangular in the basis (Cj) of Lie a i g (A sub ; k). □ 


Remark 3.8. One may refine this proof to predict the eigenvalues of T. Let 
7 i(x),..., Juj(x) £ k be the eigenvalues of i ag (x). The above reasoning shows the ex- 



/ Lfix) 


0 


istence of P\ £ GL n (C), such that P\Afc a ,fix)P 1 1 =: 




, where for 


\ 0 


Lu(x) y 



1 ^ i ^ u), Lfix) is a matrix in coefficients in k, with only one eigenvalue 'Ji(x). 


In the proof of Lemma \3. 71 we have proved the existence of a basis of constant matrices, 
such that the matrix of the linear map 4^ in this basis is block-diagonal, with blocks that 
are upper-triangular matrices corresponding to convenient restriction of the linear maps 
'kjj : Xi j >—> XijLi(x) — Lj(x)Xij. For 1 ^ i,j < w, the map 4/jj admits only one 
eigenvalue that is equal to 'yfix) — 7 fix) £ k. Then, the eigenvalues of 4' are of the form 
{li( x ) ~ lj(x), 1 ^ i,j w}. Now the diagonal blocks are symmetric powers of Ai rec i(x); 
the latter has an abelian associated Lie algebra and is triangular. It follows that the 'yfix) 
are linear combinations (with integer coefficients) of the eigenvalues of Ai re fix), so that the 
eigenvalues of *k also are linear combinations (with integer coefficients) of the eigenvalues 
of A\ re( i(x). 

3.4. Decreasing the Dimension of Lie(A;k). We refer to 112.31 113.11 and 113.21 for 
the notations and definitions used in this subsection. The aim of this section is to 
find a gauge transformation P(x) such that Y'(x) = P(x)\A(x)\Y(x) is in reduced 
form. Thanks to Corollary 13.51 it is sufficient to compute a gauge transformation 
P(x) £ jld n + B(x), B(x) £ Lie a i g (A sub ; k) j such that, for every gauge transformation 

Q(x) £ jld n + B(x),B(x) £ Lie a ig (A sub ; k) j, we have Lie(P[A];k) c Lie ^Q[P[A]];kj. 

The k-linear adjoint map <k = [•,A dia9 ] : Lie a i g (A sub ; k) -> Lie a i g (A sub ;k) has its 
eigenvalues Ai(a:),..., X K (x) in k (see Lemma [3.71) and its minimal polynomial has the form 

K, 

n#pO = n ( x - A i(x)) mi , with rrn £ N*. 

i= 1 

For each eigenvalue A fix), we let E\. := ker ((\k — Aj(a;)Id (J ) m ') be the corresponding char¬ 
acteristic space. So we have the standard decomposition Lie a i g (A sub ; k) = @1 = ^E \.. Of 
course, the E\. are 4Mnvariant subspaces. Now Lie a i g (A sub ; k) is also a 4Mnvariant sub¬ 
space of Lie a i g (A sub ; k). As the E\. have each a basis formed of constant matrices (Lemma 
13.71) . Proposition 13.61 implies that we thus have 

K 

Lie a i g (A sub ; k) = ((f) ( Ll\. Lie a j g (A sub ; k)J . 

i= 1 

In the reduction process, we may (and will) hence perform a reduction on each E\ t 
separately. So, without loss of generality, we now assume that \k has one eigenvalue 








16 


EFFECTIVE MORALES-RAMIS-SIMO THEOREM 


A(x) e k and ILj(AI) = ( X — A(x)) m , for some raeN*. 

As above, we let E\ := ker (('ll — A(x)Ido-) m ) and, for i e {0 let : = 

ker ^(T — A(x)Id (J ) t j. We have the standard flag decomposition E\ = 3™i E^/E^ 3 
And, last, we recall that for M(x) £ E^/E^ X we have 


(3.3) T(M(x)) = A(x)M(x) + M(x), with M(x) e e\'{ 1} . 

3.4.1. Reduction on One Level of a Characteristic Space. Let us first pretend that we know 
a basis C \,... ,Ct of E^'^/E^' 1 ^ (formed of constant matrices Ci, this is possible due to 
lemma [3~71) such that C s+ i ,... ,Ct form a basis of g(k) n (e^/E^' 1 1 M . This means that 
Ci,... ,C S could be “removed” by a gauge transformation. 

t 

We decompose A(x) as A(x) = A(x) + ^ Gq(x)Ci, where A(x) £ E^" 3 Our gauge 

2 = 1 

t 

transformation matrix is of the form P{x ) = Id n + X /j(x)Cj with fi(x) £ k. As ’L(Ci) = 

2 = 1 

A (x)Ci + Ci, with Ci £ E^ n 1 \ we apply Proposition 13.61 to obtain: 

t t 

P[A] = A(x) + X fi{x)Ci + X ( Mx) + X(x)fi(x) - /-(x)) Q. 

2=1 2=1 

We see that, in order to achieve reduction in E^/E^™ we should have 

/•(x) = A(x)/i(x) + aj(x) for all i e {1,... , s}. 

In other words, the differential equation y\x ) = A(x)y(x) + cq(x) should have a rational 
solution for each i e { 1 ,... , s}. 


In practice, we do not know the Ci nor the a,(x) so we now show how to compute them. 
Let B\,... ,E>t denote a basis of E^/E^ 1 ^, formed of constant matrices. We will find 

candidates for the Ci by computing which combinations of the E>i may be “removed” from 

t 

A{x) by a gauge transformation as above. We decompose A(x) as A{x) = A(x) + ^ bi(x)Bi. 

2 = 1 


t 

There exist (yet unknown) constants c t j such that Bi = ^ CijCj, so that: 

3 =i 


A(x) 


t 

A{x) + X bi(x ) 
2 = 1 



/ t 


A{x) + X X 


C 


3 - 


3 =1 V=1 


So, the calculation from the previous paragraph shows that there should exist (jj[x) £ k 

t 

such that, for j e {1,..., s}, g'j(x) = X(x)gj(x) + X Cijbi(x). The way to find s, the gj(x) 

2 = 1 

and the Cjj is given by Lemma [3.91 (which is proved here for convenience but is well known 
to specialists). 
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Lemma 3.9. Let X(x),bi(x),... ,bt(x) be dements of k. The set of tuples (g(x),c\,... ,ct) € 

t 

k x C such that g'(x) = X{x)g{x) + ^ Cibi(x) is a C-vector space. Moreover, one can 

i= 1 

effectively compute a basis of this vector space. 


Proof. Let Lb be the linear differential operator of order t whose solution space is spanned 
by b\(x ),..., bt(x). Let L := Lb • (^ — A(x)), where the product is the composition, 
i.e. the usual product in the non-commutative Ore ring k[^]. One readily sees that a 
function g(x) e k satisfies L(g(x)) = 0 if and only if L^g'(x) — A (x)g(x)) = 0, i.e. if there 

~t 

exist constants Cj £ C such that g'{x) — A (x)g(x) = ^Cjfrj(x). Hence, the set of tuples 

2 = 1 
t 

(, g(x ), ci,..., ct) e k x C such that g'{x) = A (x)g(x) + ^ Cj6j(x) is isomorphic with the set 

2=1 

of rational solutions g{x) of L. The latter is a vector space whose basis can be effectively 
computed, see H2.3.1I □ 


Lemma 13.91 allows us to compute easily, see H2.3.11 a dimension s 6 N and a basis 




j=l..s 


of elements in k x C such that the equation y'{x) = A (x)y(x) + 


y t Cijbi(x) has a rational solution y(x) = gj{x). The unknown functions ai(x) that we 
2=1 

t 


were looking for are thus given by a,(x) = ^ c l , J b l (x). 

2 = 1 

Via the incomplete basis theorem, we construct a constant invertible matrix Q e GL t (C) 
whose first s columns are the We may view Q as the base change matrix from the 

basis (Hj)* =1 of /E^ n ^ to a new basis (Cj)*- =1 of E^/E^ n Let jij denote the 

entries of Q 1 . 


Lemma 3.10. Let seN, (gj(x))j = i s > and ffi.j) be computed as in the above paragraph. 

S t 

For i 6 {1,... ,t}, let f(x) := ^ 7 ij9j{ x )- Finally, let := Id n + fi(x)Bi. Then 

j =1 i=l 

is a partial reduction matrix, in the sense that 

(3.4) Lie a i g (pW [A]-, k) n (e^/E^) = 0 (k) n (e^/E^) . 

t 

Furthermore, for all Q(x) := Id n + ^ hi(x)Ci with h s+ i(x),..., h t (x) 6 k, we have 

i=s+ 1 


Lie(pJ m) [H];k) = Lie [A]]; k) . 

Proof. We apply the first point of Proposition 12.11 (because G is connected, see the proof 
of Lemma 13.41) to deduce that g(k) c Lie a i g [A\: kj . Then, 

s(k) n (Er ) /E t p~ 1> ) s Lie^ (pf>M;k) n (e^/E^) ■ 


(3.5) 
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We want to prove the equality. By construction, C\,... ,C S vanish in the construction 
of [A] so that C s+ 1 ,... ,Ct now form a basis of Lie a i g [A]; n (^E^ n ' > / E^ n . 

t 

Due to Theorem 13.31 there exists R{x ) = Id n + ^ hi(x)Ci + R(x), with hi(x) £ k, 

z=s+l 

R(x) e E ^, such that 

(3.6) fl (k) n (Pf V4 m_1) ) = Li eaig (R^ [A]]; k) n • 

But by construction, we have the inclusion 

(3.7) Lie(Pi m) [A]; k) n (e^/E^) ^ Lie (P[P| m) [A]]; k) n (pf } /E^) • 
Combining (13.51) . (13.611 and (13.711 proves (13.4H . 

t 

Let Q(x) := Id n + ^ hi(x)Ci with h s +i(x ),..., ht(x) £ k. By construction, we have 

2=S+1 

(3.8) Lie(p| m) [A]; k) n (e^/E^) = Lie (q^ [A]]; k) n (e^/E^) ■ 

Let Cj := d'(C'j) — A (x)Cj. We use (j3.3|) and the fact that 'L is k-linear plus Proposition 
13.61 to deduce the existence of A(x) £ Lie(p{"^ [A]; k) n ^E^ n ' 1 /E^ n ^ such that 

(3.9) P^ m \x)[A(x)] ~ Q(x)[P^ m \x)[A(x)]] = A(x) + hi(x)Ci. 

i=s +1 

Let j e {s + We know that Cj £ Lie(p| m [A]; k). By definition, the matrix 

Cj = ^(Cj) — A (x)Cj belongs to Lie(P|”^ [A]; k) n E^" 1 ■*. Due to (13.8p . it also belongs to 

t 

Lie (q\P^ [A]]; kj n E^ n l> . Then, ^ hi(x)Cj belongs to Lie(P^ m i[A]; k) n E^ n L 

2=S+1 

and Lie (QlP^ [A]]; k^ n E^ n 1 \ We combine this fact and (13.91) to deduce 

Lie(p| m) [A];k) n = Lie (Q[P| m) [A]]; k) n Pf“ 1} . 

Combining (j3.8[) and this equality, we find the result. □ 

3.4.2. The Full Reduction Procedure. The reduction procedure now is easy to establish by 
iterating the above process. By assumption, all variational equations of lower order are in 
reduced form and have an abelian associated Lie algebra. 

Choose an eigenvalue A(x) £ Spec(T) of the adjoint map T = [•,A ( fj as ]. Let E\ := E^ 
be the corresponding characteristic space. Let l := m. 

Compute a constant basis of E^ / E^ L an d compute the partial reduction ma- 

t 

trix P® := Id„ + S fi(x)Bi as in Lemma 13.101 Perform the transformation A{x) : = 
2=1 

P\ ( x ) [A(a:)], let l := l — 1 and iterate this paragraph until l = 0. 
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m 

When all these successive steps are performed, let P\(x) := J (x). Note that, by con- 

1=1 

struction, the matrices P^ ( x ) commute pairwise so the order does not matter in the latter 
product. 

Now perform this for all eigenvalues A(x) £ Spec(T). The resulting matrix is a reduced 
form. 


Theorem 3.11. Using the algorithm and notations of the above paragraph, let 
P(x) : = ]~~[ P\(x) and A Pired (x) := P(x)[A(x)\. 

X(x)e Spec(^) 

Then the system Y'{x) = A Ptred (x)Y(x) is in reduced form and P(x) is the corresponding 
reduction matrix. 


Proof. Define A su b(x) as the off-diagonal part of A Pi red(aO as in the rest of this section. 


Pick any matrix H{x) e Lie a i g (A su b;k 


n 


(0 / i) 


E \t E . 


A 


for some A(x) £ Spec('L), for 


some integer l. Let Q(x) : = Id n + H{x). Then, Lemma 13.101 implies that we have the 
equality Lie(A P;red ; k) = Lie ^Qp4p,red]; kj. Now, Lemmas 13.11 and 13.21 show that any 

matrix in jld n + B(x),B(x) e Lie a i g ^A su b; kj | is a product of matrices Id n + H(x) 
of the above form. It follows that, for every gauge transformation Q(x) in the set 
Id n + B(x),B(x) e Lie a i g C4 sub ;k) j, we have Lie(A Pired ; k) = Lie (^Q[T Pired ]; kj . So, 

Corollary 13.51 shows that the system Y'{x) = d. p red (x')T (x) is in reduced form and P(x) is 
the corresponding reduction matrix. □ 


4. Back to the Morales-Ramis-Simo Integrability Criterion 

4.1. Reducing the First Variational Equation. Initially we assumed that the first 
variational equation had been put into reduced form and had an abelian associated Lie 
algebra. However, the procedure described in this paper can be also used to put the first 
variational equation into reduced form, i.e. to apply effectively the original Morales-Ramis 
integrability criterion. This allows us to recover the reduction method established by two 
of the authors in [AMW12| . 

First, factor the first variational equation, i.e. compute an equivalent lower block- 
triangular form differential system, (see e.g. |PS03| ). Then, apply a reduction procedure 
to the irreducible blocks on the diagonal (for example the one of Aparicio-Compoint-Weil 
from |AMCW13j ). This will put these blocks in diagonal form (maybe after an algebraic 
extension); otherwise we have an obstruction to integrability (Boucher-Weil criterion, see 
[BW031 IMRR10| ). If the blocks have dimension 1 or 2, then a faster method using a variant 
of the Kovacic algorithm is given in |AMW12| . 

Once this is done, the method of this paper allows us to reduce the lower triangular 
blocks, thus putting the first variational equation into reduced form. 

4.2. The Effective Morales-Ramis-Simo Integrability Criterion. The Morales- 
Ramis-Simo integrability criterion states that if one of the variational equations of a Hamil¬ 
tonian system has a differential Galois group whose Lie algebra is not abelian, then it is 
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not (meromorphically) Liouville integrable. For p e N*, let Y'(x) = A p (x)Y(x) be the vari¬ 
ational equation of order p, let G p be the differential Galois group of Y'[x) = A p (x)Y(x) 
and let g p be the Lie algebra of G p . 

As we have seen in 114.11 we may use the procedure of f|3] to put the first variational 
equation Y'{x) = Ai(x)Y(x) in reduced form. If gi is not abelian, which can be checked 
easily, then the original Morales-Ramis integrability criterion fails. Let p ^ 2 and assume 
that, for all me {1,... ,p — 1}, we know a gauge transformation matrix P m (x ) such that 
P m (x)[A rn (x)] is in reduced form, i.e. Lie a i g (P m [A m ]) = g m . We further assume that each 
g m is abelian. Then, see H2.21 the p th variational equation is of the form 


sym p (Ai(x)) 

0 

S p (x) 

A P - i(x) 


consider (see 


Let Q(x) := | 

Sym p (Pi(x)) 

0 

0 

Pp-i(x) 

sym p {A^ red {x)) 

° ^ 


S(x) 

A-p— l 5 red(^) J 



and 


Let P(x) be the gauge transformation that we have computed in 113.41 Then, 


A Pi red(x) ■= P{x)[A(x)] = P(x)[Q(x)[A p (x)]] 

is in reduced form. If g p is not abelian, which now can be easily checked, the Morales- 
Ramis-Simo integrability criterion fails. If g p is abelian, we may iterate the same procedure 
in order to put Y'{x) = A p+ \{x)Y{x) in reduced form. 

To summarize, for any p^ 2, we are able to put the successive variational equations 


Y\x) = Ai(x)Y(x ),..., Y\x) = A p (x)Y(x) 


in reduced form or prove that one of the gj is not abelian. 


4.3. A Simplified Reduction Procedure. In view of the applications of this reduction 
procedure to the Morales-Ramis-Simo integrability criterion, we have the following shortcut. 
We refer to 1}2] and 1J3] for the notations used in this subsection. The Morales-Ramis-Simo 
integrability criterion implies that, if the Hamiltonian system is integrable, once our reduced 
form from Theorem 13. Ill is computed, g p should be abelian for all p e N*. With Lemma [3.11 
we find that this is equivalent to saying that the resulting adjoint map Tred = [•,A c n ag ] 
should be the zero map (because Lie a i g (A su b) is always abelian and Lie a i g (Adi ag ) is assumed 
to be abelian). 

So, when performing the reduction, any characteristic space E\ corresponding to a non-zero 
eigenvalue A(x) 6 Spec('k) must vanish. Also, for A = 0, all (for l > 2) must vanish too. 
As a consequence, if one is only interested in finding an obstruction to integrability but not 
necessarily a reduced form, the reduction step in H3.4.1I can be significantly simplified. 

Indeed (we use the notations from H3.4.1I) . instead of the equation with parametrized 
right-hand side in Lemma [3.91 it is enough to look for a rational solution gi(x) to each equa¬ 
tion y'(x) = A (x)y(x) + bi(x). If any of these equations does not have a rational solution, 
then the adjoint map \R rec i of the reduced form will still have the non-zero eigenvalue A(x), 
hence yielding an obstruction to abelianity of the associated Lie algebra. 

Otherwise, the partial reduction matrix of Lemma 13.101 is easier to compute: just let 
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t 

P^\x) : = Id n + gi(x)Bi, compute (x)[A(:e)], compute a basis ( Bi ) of the new 
i =1 

space E^ n and iterate this reduction as in £13.4.11 Do this for all non-zero eigenvalues of 
T. For the zero eigenvalue, proceed similarly for the (for all l > 2). Note that since 
A = 0, the problem is slightly easier. Indeed, using the notations from H3.4.11 we only have 
to check whether every admits a primitive gi(x) 6 k. If any of the b,(x) does not admit 

a primitive in k, we obtain an obstruction to abelianity. Otherwise, the partial reduction 

t 

matrix will be P^\x) : = Id n + ^ gi(x)Bi. If at this stage the process has not stopped, the 

1=1 

partially reduced matrix has an associated Lie algebra which is abelian so the application 
of the Morales-Ramis-Simo integrability criterion now makes it necessary to go to the next 
higher variational equation. 

We may even iterate the process to the next variational equation without finishing the re¬ 
duction: the only assumption that was used in our algorithmic construction was that the 
Lie algebra associated to the previous variational equation was abelian. However, this is 
not very satisfying and one should, at this last step, compute the reduced form by applying 
Lemma 13.91 until the final case A = 0 and m = 1. Since A = 0, the computations here are 
slightly easier. 


5. An Example 

Consider the Hamiltonian with potential given by El := ^Pi + ^P 2 + V, where the potential 
is given by 

V = ( 9 ?i + ( li) ■ 

% 

This potential appears at the end of the |CDMP10| where the authors present it as a case 
where their necessary conditions are all satisfied so that one could guess that the system 
might be integrable, but it is left as an open case. The corresponding Hamiltonian system 
is 

( x h) ■ < 


Using the method of Darboux points and homothetic solutions (see |CDMP10| or the papers 
by the same authors or the papers by Combot in the references), we find a (rather obvious) 
pencil of particular solutions 

qi = A x,q 2 = \/"-3</i = \/^3Ax,pi = q\ = A ,p 2 = V~Spi = \/-3A. 

To simplify the expression of later results, we choose A = 4 i3 3 / 4 m £_ 1 ; the pencil now 
depends on a free parameter m and our particular solutions are: 


Qi 

<h 

Pi 

P2 


Pi 

P2 

o <?2 (3<?l 2 +<?2 2 ) 

O-t- 

Ql 

3 3gi 2 +q2 2 


4z3 3 / 4 



12-3 1/4 



4 z3 3 / 4 




i 

m 2 — 1 ’ 


P 2 


12-3 1/4 




i 


qi = 


m? — 1 
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5.1. First Variational Equation. The first variational equation has matrix A\ with 

/ 0 0 ! 0 \ 


A i = 


V 


0 

3 m 2 — 1 

8 3? 

-«V3(m 2 -l) 

8x' 2 


0 

-i-\/3(™ 2 —l) 
83 ? 

_ m 2 — 1 
8x 2 


0 1 
0 0 
0 0 / 


Following H4. 11 we find a reduction matrix for this first variational equation: 

/ x 1 i\J 3 1 \ 

—ixy/3 —i\J 3 1 —*/3-\/3 

^ q i/2-\A(rn-{- 1) _1/2 m ~— 

o 


Pi = 


V -zV3 0 1/2 2*1 

So, the reduced form of the first variational equation is 

/ 0 0 0 0 

0 0 0 0 


i/6\A(m— 1) 




A 


1 .red 


o o 


1 m +1 


2 a: 

\ 0 0 0 


0 

1 m— 1 


/ 


The associated Lie algebra is one dimensional (and abelian). We hence turn to the second 
variational equation. 

5.2. Second Variational Equation. The matrix of the second variational equation is 

A 2 (x) = 


sym 2 (A!(x)) 

0 

S 2 (x) 

A i(x) 


We start with the partial reduction matrix 

Q2,i{x) = 

to obtain 


Sym 2 (Pi(x)) 

0 

0 

Pi(x) 


A(x) := Q 2 ,i(x)[A 2 {x)] = 


sym 2 (A ljred (x)) 

0 

S 2 ,i(x) 

^-1 ,red(*^) 


where 


#2,1 (a;) = c 2 


/ 0 0 
0 0 
0 0 
\ 0 0 


0 

0 

l 

xm 

1 

mx 2 


o 

0 

j 

m 

1 

xm 


o 

0 

0 

0 


0 

0 

1 

mx 2 

_1 

mx 3 


o 

0 

1 

xm 

1 

mx 2 


1 

x* 


10/3 iV 3 10/3 iV3 


-1 

10/3 iV 3 


-10/3iV3 -10/3 iV3 -10/3 W 3 


1 


with c 2 = — (1 + i) (m 2 — 1) 3 '“ V2 -34. 
48 v ’ 


The off-diagonal Lie algebra hswfo is generated by four matrices and calculation shows that 
it has dimension 10. The matrix T of the adjoint action |/4diag, # ] on f ) su i, has eigenvalues 
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[—3/2, — 2 / 2 , —/ 2 ,0, / 2 , 2 / 2 , 3 / 2 ], where f 2 := (eigenvalue of ^4i jT . erf ), and is diagonaliz- 
able. Applying our algorithm produces the reduced form 


A 


2 ,red 


where 


S 2) red{x) C 2 


(*) y 

~ j 

&2,red 

(x) 


^4-1 ,red(^) 

( ° 

0 

0 

0 

0 

000 

0 

0 

0 

0 

0 

000 

0 

0 

1 

xm 

0 

0 

000 

V 0 

0 

0 - 

1 

xm 

0 

000 


0 0 \ 

-± 0 

X 

0 0 
0 0 / 

is a constant matrix) 


This follows from the 


We remark that Lie(A2, re d) is one-dimensional (because xA 2t red 
whereas its algebraic envelope Lie a i g (A 2:re d) is two dimensional, 
fact that an algebraic Lie algebra contains both the semi-simple and nilpotent part of each 
of its elements. This can also be seen by solving the reduced system. This is now very 
easy and the Picard-Vessiot extension is k(x — ~ ; ln(x)), which has transcendence degree 
two over k. A simple calculation (or a look at the Picard-Vessiot extension) shows that 
Lie (A 2 r ed) is again abelian so we may proceed to the third variational equation. 


5.3. Third Variational Equation. We do what we did for the second variational equation; 
we perform the first partial reduction on the diagonal to obtain 


A(x) := Q 3 i i(x)[A 3 (x)] 


sym 3 (A ljred (x)) 

0 

£ 3,1 (*) 

^4-3 ,red(*^) 


The off-diagonal Lie algebra \) su b now has dimension 33. The matrix 'L of the adjoint 
action [Adi ag ,*] on is no longer diagonalizable. Letting again f 2 := the minimal 

polynomial Ll^(V) of T is 


A 2 (X - f 2 ) 2 (A + f 2 f (A - 2 f 2 f (A + 2 f 2 f (A - 3/ 2 ) (A + 3 f 2 ) (A - 4/ 2 ) (A + 4 f 2 ) . 


Our reduction procedure turns £ 3,1 (x) into <S 3 >re< i(x) := —M 3 , where M 3 is the matrix 


( 0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
V 0 0 


0 

0 

2_ 

m 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

__2 

77 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

J_ 

m 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

_ 2 _ 

m 

0 

0 


0 

-1 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 


0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 

0000 


0 

0 

0 

0 

-1 

0 

0 

0 

0 

0 


0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 


0 

0 

0 

0 

0 

-2 

0 

0 

0 

0 


0 \ 

0 
0 
0 
0 
0 
0 
0 
0 
0 / 


The associated Lie algebra Lie a ig{A^ tre d) is still two-dimensional and is still abelian. Ac¬ 
tually, the reduced system is easily solved and its Picard-Vessiot extension is the same as 
that of VE 2 so they still have the same (abelian) differential Galois group. 
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6. Conclusion 

The reduction procedure established in this paper gives an effective version of the 
Morales-Ramis-Simo criterion in the sense that it allows us to effectively test whether 
an p-th variational equation has an abelian Lie algebra. However, when the first p — 1 
variational equations have an abelian Lie algebra but the p-th. does not, there is no known 
way to measure a priori which p would be needed. So, one may apply the reduction 
iteratively to higher and higher order but there is no criterion for determining when to 
stop. Also, when all variational equations have an abelian Lie algebra, the system could 
still be non-integrable (but one would see this on the variational equations along another 
particular solution). 

This reduction procedure will also allow further study of how the dimensions of the Galois 
groups of the successive variational equations evolve, both in integrable and non-integrable 
situations. 


The reduced form may also be combined with the methods of [AMBSWlll ISiml4b| 
for finding Taylor expansions of first integrals. Once the system is in reduced form, 
the results of |AMCW13| show that the Taylor expansions of a first integral, along the 
particular solution T, have constant coefficients. So, once the system is in reduced form, 
the (eventual) first integrals are easily found. In that sense, our reduced forms appear as 
pre-normal forms along L. Pushing the reduction further to develop a normal form theory 
would be a natural development. 


The concepts of variational equations are mutatis mutandis the same for general 
(non-Hamiltonian) dynamical systems (see e.g. |Cas09| or (CW15| . where several notions 
of variational equations are compared). The notion of Liouville integrability may be 
generalized to these contexts by Bogoyavlenskij integrability: the notion of involution 
of first integrals is replaced by the (equivalent) notion of commuting vector fields, 
see |AZ101 IBC051 Bog96 ,CB97| . The Morales-Ramis-Simo theory is generalized in 
f | AZ1 Ol ICas09| l to any kind of ordinary differential systems. The reader may remark 
that, in this paper, we essentially never use the symplectic structure of the Hamiltonian 
system from which we started. Hence, the reduction methods that we developed in 
the (symplectic) Morales-Ramis-Simo context extends naturally to any Bogoyavlenskij 
integrable differential system. 


Our reduction procedure is interesting in its own right because it applies to other kinds of 
"solvable" situations that can be found in the context of differential Galois theories. Indeed, 
consider a differential system of the form Y' = A(x)Y where A(x) has the form 


A(x) = 


Assume that the block-diagonal part 


Ai(x) 

0 

S(x) 

M{x) 

Ai(x) 

0 

0 

A 2 (x) 


is in reduced form and has an 


abelian associated Lie algebra. Our reduction procedure readily extends to this (slightly 
more general) situation and puts the system into reduced form. In particular, it may be 
viewed as a way to pre-simplify the solutions. 
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Last, we mention the case of diagonals with a non-abelian Lie algebra. In |CW15j . Casale 
and Weil develop a similar reduction technique to a family of systems in the above form but 


where 


A^x) 

0 

0 

A 2 (x) 


has a non-abelian Lie algebra. Mixing these ideas and the ones 


developed in this work may provide a way toward a reduction method for general reducible 
linear differential systems. 
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